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Abstract 

^ . The properties of isospin-asymmetric nuclear matter have been investigated in the frame- 

Ch ■ work of the extended Brueckner-Hartree-Fock approximation at zero temperature. Self- 

consistent calculations using the Argonne V14 interaction are reported for several asymmetry 
^ I parameters /3 = ^^-^ ranging from symmetric nuclear matter to pure neutron matter. The 

5^ \ binding energy per nucleon fulfills the law in the whole asymmetry range. The symmetry 

energy is calculated for different densities and discussed in comparison with other predic- 
tions. At the saturation point it is in fairly good agreement with the empirical value. The 
present approximation, based on the Landau definition of quasiparticle energy, is investigated 
in terms of the Hugenholtz-Van Hove theorem, which is proved to be fulfilled with a good 
accuracy at various asymmetries. The isospin dependence of the single-particle properties is 
discussed, including mean field, effective mass, and mean free path of neutrons and protons. 
The isospin effects in nuclear physics and nuclear astrophysics are briefly discussed. 
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1 Introduction 



Within the general interest for the equation of state ( EOS ) of nuclear matter in nuclear 
physics as well as in nuclear astrophysics, increasing attention is currently paid to the isospin 
degree of freedom. 

The EOS of isospin asymmetric nuclear matter plays a central role for our understanding 
of astrophysical phenomena like supernova explosions, neutron stars structure. X-ray bursts, 
neutron stars merging and possibly 7-ray bursts. The study of asymmetric nuclear matter 
represents also the first step for a microscopic theory of the structure of nuclei far from the 
valley of beta stability. This "terra incognita" is going to be explored in the near future 
thanks to a new generation of experimental facilities with high intensity radioactive ion 
beams. Moreover, dynamical simulations of collisions between neutron-rich nuclei show that 
the main reaction mechanisms including fragmentation are quite sensitive to the density 
dependence of the nuclear symmetry energy ^. Such calculations mainly make use of 
phenomenological Skyrme-like forces where the symmetry energy at high density can also 
be in strong disagreement with the one extracted from the microscopic predictions. 

On a microscopic basis the EOS of asymmetric nuclear matter has been studied within 
the variational approach 0, ^, Q as well as relativistic [|, |^, H, ^ [10|, |ll| and non-relativistic 



121 Brueckner-Bethe-Goldstone (BBG) theory. Within the Brueckner-Hartree-Fock (BHF) 



approximation to the BBG theory a systematic study of isospin effects on the EOS of asym- 
metric nuclear matter has been carried out in Ref . , where a separable version of the 
Paris potential Jl^ was adopted to describe the two-body nuclear force. 

Beside the bulk properties (EOS), the authors of Ref. []T^ focused also on the single- 
particle (s.p.) properties of neutrons and protons in isospin-asymmetric nuclear medium. 
The neutron and proton s.p. potentials were calculated [|l^ to the lowest order in the 
Brueckner reaction matrix (BHF approximation), using the so-called continuous choice |]T3|. 

Motivated by the renewed interest in this subject, in the present paper, we report an 
extension of the calculations of Ref. |jl2[ along the following lines. First, in the calculations 
we make use of a different realistic nucleon-nucleon ( A^A^ ) potential, i.e., the full Argonne 



Vi4 potential |16], which enables us to take into account a larger number of partial waves 



with respect to the calculation [T^ with the separable Paris potential. These additional 
partial waves ( 3 < L < 6 ) give a non-negligible contribution both to the EOS and the 
nucleon mean field, expecially in the high-density region, which is relevant for applications 
in astrophysics as well as heavy ion physics. 

Second, the Bethe-Goldstone equation is now solved for the complex G-matrix. This 
enables us to calculate the complex nuclear mean field and some closely related quantities 
such as the optical potential and the mean free path. 

Third, according to the Landau definition of quasi-particle energy (for an extended dis- 



1 



cussion see Ref. [|I^) in the calculations of the mass operator (nucleon self-energy) and 
single-particle properties, we go beyond the BHF approximation by including some higher- 
order correlation contributions. In particular, in the present work, we include the so-called 
rearrangement term M2 which is a second order diagram in the G-matrix and accounts for 
particle-hole excitations in nuclear matter ground state. Next we consider also the renor- 
malization contributions of the third and forth order in the G-matrix, which account for the 
partial depletion of the neutron and proton Fermi seas due to the nuclear correlations p 



It has been shown, in the case of pure neutron matter 19] and also symmetric nuclear mat- 



ter that the new terms give a large contribution to s.p. properties like the mean field 
and the nucleon effective mass. We will refer to the present approach to compute nuclear 



s.p. properties as extended Brueckner-Hartree-Fock ( EBHF ) approximation |T^, ^0 . 

As is well known, the BHF approximation largely violates the Hugenholtz-Van Hove 
( HVH ) theorem |]2T|, which basically measures the consistency of a given order of ap- 
proximation in a perturbative approach. In symmetric nuclear matter, the inclusion of the 



rearrangement contribution greatly improves the fulfillment of the HVH theorem In 
the present paper, we study this problem in the case of asymmetric nuclear matter within 
the EBHF approximation. 

2 EBHF and Nucleon Self-energy for Asymmetric Nu- 
clear Matter 

In this section the formalism of the Brueckner-Bethe-Goldstone ( BBG ) theory is described 



for the case of asymmetry nuclear matter |12, 23]. The proton and the neutron Fermi 



momenta are related to their corresponding densities pp and pn through the relations 

k'-F = [—{I + m'^', 

where p = pp + pn is the total density, and /5 = (pn — Pp)/p the asymmetry parameter 
determining the neutron excess (from now on we assume pn > Pp). 

The starting point in BBG theory, is the Brueckner reaction matrix G, which in the case 
of asymmetric nuclear matter depends also on the isospin components of the two colliding 
nucleons. The G-matrix satisfies the Bethe-Goldstone equation, 
^/ ^ N , ^ \kik2)Q{ki,k2){kik2\ 

l^,^-<ki)-e{k2)+tv 

where vnn is the two-body nuclear interaction and uj the starting energy. Here k = {k, a, r) 
denotes s.p. momentum, z-components of spin and isospin, respectively. 
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The G-matrix can be considered as an in-medium effective interaction between two nucle- 
ons. The surrounding nucleons renormahze the bare NN interaction via the Pauh bloclcing 
and the nuclear mean field. The Pauli operator, defined as 

Q{h,k2) = [l-n{h)][l-n{h)] , (2) 

prevents two nucleons in intermediate states from scattering into states inside their respective 
Fermi seas. By n{k) we denote the Fermi distribution function, which at zero temperature 
is given by the step function 9{k — kp) ( uncorrelated ground state ). The s. p. energy [] 

eik) = ^ + U{k) , (3) 

appearing in the energy denominator of Eq. (|l]), involves the auxiliary potential U{k), which 
controls the convergence rate of the hole-line expansion. Within the BHF approximation the 
neutron and proton s.p. auxiliary potentials are calculated from the real part of the on-shell 
antisymmetrized G-matrix, via the relation 

U{k) =J2n{k')Re{kk'\G{e{k) + e{k'))\kk')A ■ (4) 
k' 



Here we adopt the continuous choice |T5| for the auxiliary s.p. potential. In this context, 
the auxiliary potential has the physical meaning of the mean field that each nucleon feels 
during its propagation between two successive scatterings. 

In the BHF approximation, Eqs. (|I]),(|^), and are solved self-consistently for given 
total density p and asymmetry f3. Then the energy per particle is evaluated at the lowest 
order ( two hole-line diagrams ) of the BBG hole-line expansion ( see Ref. for the case 
of asymmetric matter ) . 

2.1 Mass operator and quasi- particle energy 

One of the main purposes of the present paper is to calculate s.p. properties of neutrons 
and protons in asymmetric matter going beyond the BHF approximation. To this end we 
introduce the mass operator [ITR 



M^{k,u;) = V^{k,Lj) + iW^{k,u) , (5) 

which is a complex quantity and can be identified with the potential energy felt by a neu- 
tron (r = n) or a proton (r = p) with momentum k and energy uj in asymmetric nuclear 
matter ( hereafter, we will write out explicitly the isospin index r ). In the same spirit of the 
BBG theory, the mass operator M'^{k, u) can be expanded in a perturbation series according 
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Figure 1: Hole- line expansion of s.p. potential. 



to the number of hole lines |^ and various terms of this expansion can be represented by 
means of Goldstone diagrams a few of which are shown in Fig. 1. 

In analogy with the case of symmetric nuclear matter, the neutron and proton quasi- 
particle energies E'^{k) are the solutions of the energy-momentum relation 



27.2 



2m 



+ V^{k,E\k)) , 



(6) 



i.e., E'^{k) is obtained from the on-shell values of the real part of the mass operator. 



u 



(A;) 



+ 



n 



Figure 2: The first order hole-line expansion of neutron s.p. potential. 

To the lowest order in the hole-line expansion the mass operator is given by ( diagrams 
of Fig. 2 ) 



Ml{k,u) = E E n^' {k'){kk'\G^^' {uj + t^' {k'))\kk') A = E W'ik, 



(7) 



In this approximation the quasi-particle energy El{k) coincides with the BHF s.p. energy 
given by Eqs. (D and (|), i.e., El{k) = e^{k). 



2.2 The rearrangement contribution to the s.p. energy 

The next contribution to the perturbative expansion of the mass operator is given by the 



so-called rearrangement term M2{k,Lu) |15|. The associated Goldstone diagrams are shown 



"'^In the present work, we assume the neutron and proton rest masses equal to their average value m 
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Figure 3: The second order hole-line expansion of neutron s.p. potential. 

in Fig. 3. MJ is a second-order diagram in the G-matrix and accounts for particle-hole 
excitations in nuclear matter. Its expression, extended to asymmetric nuclear matter, reads 



k'a' 



k\k2 



where e^(fc) is the s.p. spectrum in BHF approximation, given by Eqs. and ( 
this approximation for the mass operator [ i.e., M'^{k,uj) ~ Ml{k,uj) + M2{k,uj) 
quasi-particle energy (H) is given by the approximate relation 

E^.ik) = El{k) + Zl{k)V^{k,El{k)) 



2m 



+ V{ik, Elik)) + Zl{k)V^{k, Elik)) 



where 



duo 



V,^{k,uj) + V^{k,u) 

J uj=El{k) 

is an approximation of the quasi-particle strength for asymmetric nuclear matter 

-1 



Z^{k) 



duj 



V^{k,u) 



(8) 

In 

the 



(9) 



(10) 



2.3 The renormalization contributions to the s.p. energy 

Due to many-body correlations the two Fermi seas are partially depleted, and the correlated 
momentum distributions n^{k) differ from the uncorrelated ones n'^{k) = 9{k — kp). To 
account for this physical effect, one considers the contribution M^{k,uj) ( last diagram of 
Fig. 1 ) given by [0, 

M^{k,u;) = -J2J2'^2ih'){kh'\G^-'{u; + e-'{h'))\kh')A , (12) 

t' h'cr' 
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where h' refers to "hole" state with momentum smaller than /cp, and 



- oj=e^'{h') 



is at the lowest order the depletion of neutron (proton) Fermi sea JTS], [T^, i.e., n^' {h') is the 
probability that a neutron (proton) hole-state {\h'\ < kp) is empty. Let us consider now the 



sum 



= E E [l - {kh'\G'^'{uj + e^'{h'))\kh')A 

= EE nl'{h'){kh'\G^^\u; + e^\h'))\kh')A, 



(13) 



T h'a' 



^2 i^') — ^ ~ ^2^) being the second-order approximation for the correlated momentum 
distribution. Ml{k,uj) is the so-called renormalized BHF approximation for the off-shell 
mass operator [ compare to Eq. (0) ]. 

An accurate approximation consists in using the average value of the depletion, which is 



K^' = nl'(h' = 0.75k^) . 

Then Eqs. (0) and (|T3p yield 



(14) 



(15) 



(16) 



From the similar considerations, a renormalization correction should also be brought from 
the four hole-line terms to the second-order contribution to in order to take into account 
the fact that the hole-state ki in Eq. (^) is partially empty ( see also Ref. |]TB[ for symmetric 
nuclear matter ). Along the same line of the previous correction one gets the renormalized 
M2, which is approximately given by 



MJ(A;,a;) = E[ 



MY' 



{k,uj) 



(17) 



The renormalized contributions can also be traced to the functional dependence of the G- 
matrix on the quasi-particle occupation numbers within the Landau theory of Fermi liquids. 
It can be shown, in fact, that taking the functional derivative of the binding energy (at two 
hole-line level) includes also the terms of the third and fourth order in the self-energy, the 
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effect of which has just been discussed. Taking into account all the corrections discussed 



above, from Eq. (Jg) one can get the following expression for the quasi-particle energy |T8 



2m 



+ Z^{k) ^ [-^.-'vr'ik, El{k)) + (1 - K^')Vr\k, El{k)) 



where 



vr\k,uj) + vr\Kuj) 



(19) 



u}=El(k) 



In the following we refer to this approximation for the quasi-particle energy as the extended 
Brueckner-Hartree-Fock (EBHF) approximation [O, [T^, 120 . 



2.4 Partial wave expansion and angular averaging 

After the usual angular averaging on the Pauli operator and the energy denominator PBl 
the Bethe-Goldstone equation can be expanded in partial waves, 

C^LL' (g, q\ P, UJ) = VaLL' (?, q') 



+ -E/?"WVll"(9,0 



{Q'''{q",P)) 
u-el^'{q",P)+ir] 



G-;,^,(g",g',P,a;) 



(20) 



where q = {ki — k2)/2 and P = ki + k2 are the relative momentum and total momentum, 
respectively. e'[2{q",P) = (e^(A;i) + e^' {k2)) is the angle average of the energy denominator. 
The angular-averaged Pauli operator is 

(i) for T = t' (neutron-neutron or proton-proton). 



{Q^^iq,P)) 




if > 

otherwise , 



(21) 



(ii) for T ^ t' (neutron-proton or proton-neutron), 

1, 

{Q^^'iq,P)) ' 



[min(l,^p) -hmin(l,^n)] 



if -ep,en> -1 

otherwise , 



f22) 



where 



pyi + q^-{k^y 

Pq 
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The mass operators Mi and M2 become 

1 + S.- -' — rk- 



Mr {k,Uj) 



2n 



aL 

2(l + (^r,rO 
TT^k 



^^{2J+1) [ " k"dk'smedeGl^;^iq,q,P,u + e^'{k')) , (23) 



aLL' 



^ (2 J + 1) / / qdqPdP 



[JPy2 + 2q^-k^ 



J U + e-"' (^PV2 + 2?2 _ y^2j _ grr' (^/^ _ 



(24) 



The integrations of q and P in the expression of M2 are hmited to 



?min * 



max 



max 



0, k 2(^f + ^f)) 2^^^ 



if A; < A; 



2; fvp 



V^2P - 2(A;5)2 + {kl - A;f')2, ^(A; + A;j' 



if k> kl 



k + i^{kl + kl) 



(25) 



(26) 



and 



p . 



max 



-Pmax = min 



2(A;-g), V2A;2 + 2(A;|;')2 - 4g2 

2\q-k\ 

2{k + q), {k^ + k^ 



if g<|(A; + A;j;') 
if k>l{k + k^) 



(27) 



(28) 



The angular averaging of the anti-Pauh operator EJ'^' [ki, k2) = rf {ki)rf' {k2) can be written 



as 



(i) for T = r', i.e., neutron-neutron or proton-proton 



(P--(g,P)) 



min(l, 77^) if 77,- ^ 

otherwise , 



(29) 



(ii) for T 7^ r', i.e., neutron- proton or proton-neutron 



(P-'(g,P)) 



-[min(l,77p) +min(l,?7n)] 



if ?7p > -?7n, ?7p > -1 
otherwise , 



(30) 



where rf — — 



3 Results 



We have performed a set of nuclear matter calculations for the asymmetric case within the 
EBHF approximation. Three different densities have been selected: po/2, Po and 2po, being 
Po = 0.17 fm~^ the saturation density of symmetric nuclear matter. For each density the 
whole range of asymmetry parameter ( < /3 < 1 ) has been spanned. The self-consistent 
solution of the Bethe-Goldstone equation yielding simultaneously G-matrix and auxiliary 
potential U^{k) needed five iterations to reach a satisfactory convergence. The bare potential 
adopted as input in the calculation was the Argonne V14 |jl6[ with 24 channels up to L = 6. 



3.1 Symmetry Energy 

In Fig. 4 (left panel) we report the results (symbols) for the energy per nucleon B{p,P), 



calculated self-consistently within the BHF approximation p3|. B{p,l3) is plotted as a 
function of jS"^, for three values of density. The numerical results lie on a linear fit performed 
with only the first three values of the asymmetry parameter. This proves that the empirical 
parabolic law 

B{p,f3)=B{p,0) + Esym{p)P\ (31) 

taken from the nuclear mass table can be extended up to the highest asymmetry of nuclear 
matter, in good agreement with our previous BHF calculation with the separable Paris 
potential |jl2 . 



Equation ( |3lD can be considered as the (3 expansion of the binding energy truncated 



at the lowest order. Only even powers of the asymmetry parameter f3 may occur in the 
expansion for charge-independent A^A^ interactions, such as the Argonne V14 used in the 
present work. A contribution might arise at the three hole-line order of the BBG ex- 
pansion. Unfortunately, no such calculation for B{p,/3) has been done yet. However, it has 
been shown recently that the three hole-line contribution to the binding energy of symmetric 
nuclear matter |^ is rather small within the continuous choice. Therefore, we do not expect 



a large deviation from the parabolic law after including the three hole-line contribution. A 
deviation from the parabolic law could be expected at densities higher with respect to those 
considered in the present paper p, ||, ||, |^ . 
The symmetry energy is defined as 

'd'B{p,(3) 



^sym (p) 



(32) 

J 13=0 



Due to the simple /3^-law the symmetry energy can be equivalently calculated as the dif- 
ference between the binding energy of pure neutron matter and symmetric nuclear matter: 
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Figure 4: Left panel : Total binding energy per nucleon in the range < < 1 at three 
densities as compared with the parabolic fits ( straight lines ) obtained from the first three values 
of P (0.0, 0.2, 0.4). Right panel : Density dependence of the symmetry energy of the present 
work ( solid curve ) using Argonne V14 as bare interaction in comparison with other non-relativistic 
calculations. The dashed curve is the result of the lowest order constrained variational calculation 
using Argonne V14 as bare interaction from Ref.[5]. The dotted and dot-dashed curves are the 
results of the variational approach using Argonne V14 and Argonne V14-I-UVII, respectively, taken 
from Ref.[3] 

Esym{p) = B{p,l) — B{p,0), but one would refrain from using that recipe at very high 
density. The results of our BHF calculations for Esym{p) are depicted by the continuous 
curve in the right panel of Fig. 4. In the same figure, we also show the results from the 
variational approach using the same Argonne V14 potential 0. 

The systematic disagreement displayed by the two many-body approaches has been be- 
lieved to be a shortcoming of the Brueckner approach in view of the fact that the BHF result 
lies above the variational one. However, in Ref. (and similar works), the variational ex- 
pectation value E^ar of the Hamiltonian is calculated in a diagrammatic cluster expansion 
(FHNC-SOC), which is of course truncated to some order. To estimate the convergence of 
this diagrammatic cluster expansion, we plot, in the same figure, the results of a lowest-order 
constrained variational calculation , which includes only two-body cluster contributions to 
Eyar- Moreover, the variational trial wave function used in Ref. does not contain the cor- 
relations which arise from L^, L^(ai ■ 0^2), and {L ■ SY terms of the nucleon-nucleon potential. 
Finally, spin-orbit correlations are not treated accurately, as discussed in the same paper . 
All these A^A^ correlations are included in a self-consistent way in the BHF approach. All the 
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above-mentioned approximations could give large uncertainties in the calculated expectation 
value of the energy in the high density region. The same discrepancy has also been observed 
in our previous calculations for asymmetric nuclear matter |jl2| and also in neutron matter 
calculations 

From the previous discussion we guess that the nice agreement between our calculation 
and a lowest order constrained variational calculation P] is fortuitous. On the other hand, 
the agreement up to p ~ 0.24 fm~^ with the variational calculation including three-body 
force B, also plotted in Fig. 4, is hardly understandable. 




Figure 5: Symmetry energy vs density in the BHF approximation. Left panel: symmetry energy 
(Esym), kinetic {Kgym) and potential {Vsym) contributions (only Argonne V14). Middle panel: 
potential contribution from isospin T=0 channels with Argonne 1^14 (sohd curve) and separable 
Paris (dashed curve). Right panel: potential contribution from isospin T=l channels with Argonne 
Vu (solid curve) and separable Paris (dashed curve). 

More recent versions of the NN potential do not provide any appreciable difference of the 
symmetry energy from the present calculation except the CD-Bonn potential as discussed 
in Ref. |2^. All Brueckner calculations predict the symmetry energy to increase with the 



nucleon density and no saturation is observed up to p = at variance with the 

preceeding variational results 0. In the relativistic mean-field theory this behavior is easily 
understood in terms of the p-meson exchange, which leads to a repulsive symmetry potential 
at all densities Jl], 0. In order to try to explain what happens in the nonrelativistic case, 
we report in Fig. 5 the different contributions to the symmetry energy, plotted as a function 

2 

of density. The kinetic contribution monotonically increases as pa according to the free 
Fermi-gas model. In the figure (right two panels), the isoscalar and isovector contributions 
of potential part are plotted separately. The density dependence of the symmetry energy 
is dominated at high density by the kinetic contribution, where the opposite behavior vs 
density of the potential contributions from T = and T = 1 channels results in a very flat 
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density dependence of the symmetry potential. As already found in the previous paper ||T2 



the most important contribution to the T = component is due to the deuteron ^Si-^Di 
coupled channels of the interaction, which exhibits a maximum at p ~ 0.3/m~^. This 
peak can be traced back to the behavior of the two components: the attractive ^Si channel 
dominates at low energy whereas the repulsive ^Di dominates at high energy. Two terms 
compensate each other at the energy E ^ AEp — 200MeV, where Ep is the Fermi energy 
corresponding to p 0.3/m~^. 



3.2 Single-particle energy 

For asymmetric nuclear matter the neutron mass operator is different from the proton 
mass operator M^. Moreover, as shown in Figs. 2 and 3 [see also Eqs.(7) and (8)], both of 
them can be split into two components: = M^^ + M^" for protons and M" = M"" + M"'^ 




Figure 6: Real part (upper panels) and imaginary part (lower panels) of the first-order single- 
particle potentials Mi for proton (left panel) and neutron (right panel), respectively, as a function 
of momentum for different asymmetry parameters at density p = 0.17 fm~^. 



for neutrons. In Fig. 6 the on-shell values of the real part (upper panels) of are reported 
as a function of the s.p. momentum, for different values of the asymmetry parameter (3 
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at fixed density p = 0.17fm'^. The proton mean field V[{k) = ReMf(A;) becomes more 
attractive, while the neutron mean field V{^{k) = ReM"(/c) becomes more repulsive going 
from symmetric {(3 = 0) to neutron {(3 = 1) matter. The j3 dependence of and Vf is 
almost linear and nearly symmetric with respect to their common value at /? = 0. This result 
supports from a microscopic point of view the validity of the so-called Lane potential |3C]. 

It is worth noticing that a crossing point occurs for both Vi and V{^, where the isospin 
effect on neutron and proton mean fields versus (3 is inverted. This behavior of the neutron 
and proton mean field can be understood in terms of phase-space arguments, as pointed out 
Ref . |T^ . To this end, we write the s.p. potentials and Vf in terms of their components 
V{'^' [defined according to Eq.(7)]: 

vrik) ^ 1(1 - (3)p{Gn + 1(1+ f3)p{Gn m 

V-{k) ^ 1(1 - + 1(1 + (3)p{Gn (34) 

where (G^'^) is the average value of the real part of the matrix G^*' in the proton Fermi sphere 
{\h'\ < kp), and (G^") the average value of the real part of the matrix G^" in the neutron 
Fermi sphere {\h'\ < kp). (G"") and (G"^) have similar definitions. This approximation is 
suggested by the almost linear dependence of Vi and on (3 and, in fact, is numerically 
fulfilled with a good accuracy (see also Fig. 8). The crossing point in momentum space is 
determined by the occurrence of (G^^) = (G^") for Vf and (G"^) = (G"") for at a certain 
value of the momentum which does not depend upon f3. A signature of the inversion of the 
isospin effect at the crossing point could be found in those collective observables measured 
in heavy ion collisions which are sensitive to the momentum dependence of the mean field. 

The imaginary part W[ of the mass operator M{ is due to the virtual collisions of a single 
nucleon with a neutron or a proton of the background, promoting it to a particle state. W[ 
is vanishing below the Fermi momentum kp due to the Pauli blocking. It is worth noticing 
that reducing the proton Fermi momentum implies a less Pauli blocking for protons. This 
means that high asymmetric nuclear matter is less transparent to the proton propagation. 

The second-order terms of the on-shell mass operator M2 are plotted in Fig. 7. The real 
part V2 (upper panels) gives the contribution to the mean field due to the coupling of the 
single-particle motion with the ground state particle-hole excitations. As is well known, V2 
is repulsive and reduces to a large extent the pure BHF mean field V{ which is too attractive 
compared with the phenomenological optical potential p6| . 

The imaginary part W2 plays a role complementary to W[: it describes the virtual 
collisions of a single nucleon of the background with an excited neutron or proton, making 
it to decay into a hole state. W2 is vanishing above the Fermi momentum kp. 

In order to focus on only the isospin dependence, we plot in Fig. 8 the mass operator as 
a function of /5 at = 0/m~^ except for the imaginary part of Mi for which a value of k 
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Proton s.p. potential Neutron s.p. potential 




Figure 7: Real part (upper panels) and imaginary part (lower panels) of the second-order single- 
particle potentials M2 for proton (left panel) and neutron (right panel), respectively, as a function 
of momentum for different asymmetry parameters at density p = 0.17 fm~^. 

above the Fermi momentum has to be taken. 

The first-order contribution has the hnear behavior for the real part as well as for the 
imaginary part as expected from phase-space arguments. The slope of |iyf"| is more pro- 
nounced than that of since the neutron particle-hole excitations coupled to a proton 
in a particle state are more favored than the proton particle-hole excitations coupled to a 
neutron (see also Fig. 6). 

The isospin dependence of the second-order contribution is affected by the coupling 
between the nucleon hole stats and particle-hole excitations [see the bubble in Fig.3 and 
Eq.dl)], which yields a nonlinear variation of the mixed components Mf" and M;^^ vs (3. 
The nonlinearity is much more sizeable for and PF^^, which can be easily explained as a 
phase-space effect as well, i.e., of the interplay between the neutron and proton phase-spaces 
as increasing neutron excess. 
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Figure 8: Different components of the first-order (upper panels) and the second-order (lower 
panels) single-particle potentials at density p = 0.17/m~^ for suitable values of momentum versus 
asymmetry parameter. For the real part of M\ and M2, and the Imaginary part of M2, the 
momentum is = 0, while for the real part of M\ the momentum is /c = 3/m~^. 

3.3 Fermi Energy and Hugenholtz-Van Hove theorem 

The EBHF approximation basically relies on the Landau definition of quasi-particle energy 



as showed in sec. I, whose relation to the Brueckner theory has been well established [pT |. 
Study of the HVH theorem within the EBHF approximation could provide an additional 
support to a proper definition of the quasi-particle energy and, at the same time, a more 
realistic evaluation of the Fermi energy. 

Strictly speaking, the HVH theorem concerns only symmetric nuclear matter at satura- 
tion point (P = 0), and it states that the energy per nucleon must be exactly equal to the 
Fermi energy. In the case of asymmetric nuclear matter (two-component system) at zero 
temperature, the HVH theorem can be generalized via the thermodynamic relation 



E{p, P) ^ P{p^ ^ YpEUp, 13) + Y'^El^p, (3) 



(35) 



A p 

P{p,P) being the pressure, = pp/p and F" = Pn/ p the proton and neutron fractions, re- 
spectively. The Fermi energy is calculated from the quasi-particle energy spectrum at Fermi 
surface according to Eq. (6). 
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TAB. I As a function of the asymmetry parameter (first comumn) are reported (in MeV) the physical quan- 
tities involved in the Hugenholtz-Van Hove theorem: pressure over p (second column), energy per nucleon 
(third column) and "weighted" chemical potentials of asymmetric nuclear matter in different approxima- 
tions, as discussed in the text. Y'p Z/A and F" = N /A are the proton and neutron fractions, respectively. 
The total density is p = 0.17/m~'^. 



/5 


P/P 


E/A 


P/p + E/A 


BHF 


yP£'^ + VE"^ 

BHF+M2 


yp£;^ + Y^E'^ 

EBHF 


0.0 


-5.02 


-15.92 


-20.94 


-34.27 


-28.50 


-19.28 


0.2 


-4.40 


-14.73 


-19.13 


-32.29 


-26.43 


-17.35 


0.4 


-3.27 


-11.36 


-14.63 


-26.28 


-20.74 


-12.42 


0.6 


0.08 


-5.75 


-5.67 


-16.56 


-11.44 


-4.41 


0.8 


3.76 


2.24 


6.00 


-2.40 


2.07 


7.42 


1.0 


8.91 


12.83 


21.74 


16.28 


19.67 


22.39 



In Tab. I it is numerically shown to what extent the HVH theorem is fulfilled by the 
EBHF approximation. The pressure has been calculated using the relation P{p,P) = 
p^[dEA{p, P)/dp], where Ea{p, P) = E{p,P)/A being the energy per nucleon. In the forth 
column the left-hand side of Eq. ( pS] ) is calculated for several asymmetries (density fixed at 
p = 0.17/m~^). One would notice that, despite the fact that the total density is fixed at the 
empirical saturation value, our calculated saturation point lies at higher density, because, as 
is well known, Brueckner theory with two-body force misses the empirical saturation point. 
The last three columns provide different approximations for the right-hand side of Eq. (p5|). 
The pure BHF approximation by itself is far from fulfilling the HVH theorem. Including 
the unrenormalized ground-state correlations (indicated by BHF-I-M2 in the table), where 
the Fermi energy is calculated according to Eq. (9), provides some improvement but it is 
not enough to fulfill the HVH theorem. One needs to include both the rearrangement and 
the renormalized contributions (EBHF) if a satisfactory agreement within less than 10% is 
to be attained (last column of Tab. I). This result is in keeping with the uncertainty in the 
calculation of the pressure because the binding energy curve is rather fiat as a function of 
density. 



4 Applications 
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4.1 Effective mass 



The effective mass incorporates the non-local part of the mean field which makes the local 
part less attractive for a nucleon travelling with momentum k > 0. It is defined as 



m*{k) _ k fdE^ik)' 
m m \ dk 



(36) 



The momentum dependence of m* is characterized by the wide bump inside the Fermi sphere 
due to the high probability amplitude for particle-hole excitations near the Fermi surface 



[[l5i . The effect of correlations is a flattening of the slope of the mean field around the 
Fermi energy, which implies an enhancement of the effective mass at kp with respect to 
BHF value [1^, This result is shown in Fig. 9, where an increase from 0.8 to 0.92 is 



observed for symmetric nuclear matter at the empirical saturation density. Also shown in 
the figure is the isospin dependence of the neutron (upper curve) and proton (lower curve) 
effective masses. In both the BHF and EBHF calculations, m* increases and m* decreases 
as increasing (3. Compared to the BHF approximation, the corrections of EBHF shift m* 
and m* to higher values, a feature which can be traced to the depletions of the proton and 
neutron Fermi surfaces due to the ground-state correlations. The value of m* calculated 
from EBHF approaches its BHF value as increasing j3 since the correlations become smaller. 




0.6 ' 

0.0 0.2 0.4 0.6 0.8 1.0 

(N-Z)/A 

Figure 9: Proton and neutron effective masses versus asymmetry parameter at density p = 
0.17/m-3. The solid curves are results from the pure BHF calculation, while the dashed curves 
are calculated from the EBHF (including the renormalization contributions). 
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4.2 Mean free path 



Information on the in-medium cross section or, equivalently, on the mean free path of a 
nucleon travelhng inside a nuclear medium can be obtained from the transparency of a 
nucleus measured in (e, e'p) reactions and, in general, from nucleon- induced reactions at 
low energy [^. The underlying assumption is that the behavior of a nucleon located at the 
position r'in a nucleus is the same as a nucleon in nuclear matter at density p(r). Such an 
assumption is the well known local-density approximation (LDA) The mean free path 
is intimately related to the imaginary part of the optical potential or, equivalently, to the 
imaginary part of the mean field. The latter comes from the collisions of a single nucleon 
with the background of neutrons and protons: a nucleon with momentum k > kp can collide 
with a neutron or proton of its Fermi sea and promote it to a particle state, or a nucleon 
with momentum k < kp interacting with an excited neutron or proton can make it decay 
into a hole state. The first process is related to the imaginary part of Mi, the second one 
to the imaginary part of M2, both of which have been plotted in the lower panels of Figs. 6 
and 7. But, in the case of asymmetric matter the collisions between like and unlike nucleons 
yield contributions to the mean field which are very different. 
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Figure 10: The energy dependence of proton (upper panels) and neutron (lower panels) mean free 
paths for different asymmetry parameters at three densities p = 0.085fm~'^, p = 0.17fm~^, and 
p = 0.34fm~^. Only EBHF results are reported. The different kinds of lines correspond to different 
asymmetry parameters with the same notation as in Fig. 6. 
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The mean free path A,- is given by 

h''k{E) 



2mr \lmMr{k{E),E)\ ' 

where m,- is the so-called fc-mass, and E is the single-particle energy ||3^. In Fig. 10 the 
proton (upper panels) and neutron (lower panels) mean free paths calculated within the 
EBHF approximation are shown for three values of the total density. In each panel the values 
of Xr for several asymmetries are plotted as a function of single-particle energy. The most 
relevant effect of the isospin asymmetry is the increasing deviation from the symmetric values 
(solid lines), upward for A„ and downward Ap, as increasing asymmetry. The nonvanishing 
values of neutron and proton mean free paths below their respective Fermi energies are 
effects of ground-state correlations which prevents a full occupancy of the Fermi spheres. 
Comparing with the BHF calculation it turns out that the correlation effects tend to rise the 
asymptotic value of the mean free path from about 3fm up to about 4fm at the saturation 
density of symmetric nuclear matter pO . 




Figure 11: Proton and neutron inverse mean free paths versus asymmetry parameter for three 
densities p = 0.085fm~^, p = 0.17fm~^, and p = 0.34fm~^ at a fixed single-particle energy E'^{k) = 
ISOMeV from the EBHF calculation. 

In Fig. 11 it is shown how the isospin dependence of the inverse A develops as increasing 
nuclear matter densities at a fixed value of the energy. Except for very small asymmetries 
the shift of Ap and A„ is not symmetric with respect to their common value at /5 = 0. At 
any density the slope of the neutron inverse A is less than the proton one. This effect can be 
traced to the reduction, as increasing neutron excess, of the proton particle-hole excitations 
contributing to the neutron optical potential (see Figs. 6 and 7). Moreover, the EBHF 
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Xn seems to reach the asymptotic value of pure neutron matter much faster than in the 
uncorrelated case. 

The most striking effect of the isospin-asymmetry is the sizeable reduction of the proton 
mean free path at high asymmetry. Accordingly, the nuclear surface would become more 
transparent to neutrons than protons in nucleon-induced reactions on nuclei near the neutron 
drip-line. This effect would be more pronounced at higher density as shown in Fig. 11. 

4.3 Proton fraction in /3-equilibrium matter 

The core of a neutron star is expected to be formed by an uncharged mixture of neutrons, 
protons, electrons and muons in equilibrium with respect to the weak interactions (/3-stable 
matter). The concentrations of different particles are then obtained under the requirements 



Pp = Pe + Pf, ■ (38) 

The difference between the neutron and proton chemical potentials can be expressed as 



Pn pp 



dB 



dYP 



dB 



(39) 



In the parabolic approximation, Eq. (0), for the energy per particle of asymmetric nuclear 
matter, one has 



Pn-Pp = 4E,yrn{p){l-2YP) 



(40) 



Therefore, the composition of /^-stable matter, and in particular, the proton fraction 
present at a given density, is strongly dependent on the nuclear symmetry energy. The 
proton fraction plays also a crucial role in the thermal evolution of neutron stars. In fact, if 
the proton fraction in the core of a neutron star, is above a critical value Y^^^^, the so-called 
direct Urea processes can occur If they occur, the direct Urea processes enhance the 
neutrino emission and neutron star cooling rate by a large factor compared to the standard 
cooling scenario. The critical proton fraction has been estimated to be in the range 
11 - 15%. In a recent paper [^, based on microscopic EOS of dense matter, it has been 
found that the onset of the direct Urea processes occurs at densities p > 0.54-0.65 fm~^, 
depending on the nuclear interaction used to get the EOS ( see Ref. for more details ). 

In Tab. II, we report our present calculations of the proton fraction YP{p) for /5-stable 
matter in comparison with the one obtained with the separable Paris ( see Ref. ||12[ ) and 
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variational calculation of Ref. |0 using the Argonne Vu potential plus the Urbana model 
( UVII ) three-body force. In the calculations reported in Tab. II, muons have been not 
included. 

Tab. II Proton fraction in /3-stable nuclear matter (no muons) versus the total baryonic density from 
different forces. The values reported are lO^y. The results in the second column are taken from Ref. [ p2[ . 
Those in the third column have been given by A. Fabrocini (private communication). 





Paris 


AV14+UVII 


present 


0.038 


2.75 






0.076 


2.80 


1.85 


2.40 


0.11 


3.09 


2.48 


2.74 


0.14 


3.48 


2.96 


3.03 


0.17 


3.70 


3.37 


3.32 


0.20 


4.10 


3.74 


3.50 


0.30 


4.90 


3.67 


4.07 


0.40 


5.79 


3.56 


4.57 


0.50 




3.63 


5.01 



Our purpose, in the present paper, is not an accurate determination of the proton fraction 
in dense stellar matter. Here, we aim to study how the inclusion of contributions beyond the 
BHF to the chemical potentials could alter the proton fraction in /3-stable matter. In fact, 
to solve the /3-equilibrium conditions (|37| ) and (|38|) , the shift between neutron and proton 



chemical potentials /i = //„ — yUp has to be evaluated. In Tab. Ill the neutron and proton 
chemical potentials and their difference fi, are reported for the different approximations used 
in the present work. From the results reported in Tab. Ill we see that the chemical potential, 
approximated by the Fermi energy, in the EBHF are noticeably affected by the rearrange- 
ment and renormalization contributions. However, their difference and consequently the 
proton fraction is almost unchanged with respect to the BHF approximation. The EBHF 
approximation provides neutron and proton Fermi energies, which are in better agreement 
with the empirical values extracted from the mass table of atomic nuclei |3^ than the BHF 
approximation [p!2| . 



21 



Tab. Ill Proton and neutron chemical potentials (Fermi energies) calculated in different approximations 
and compared with the symmetry energy. The results corresponding to four values of asymmetry parameter 
[3 for each of the three densities are reported. 



p 


13 




^" 


A 






A 




^" 


A 


syrn 






BHF 


BHF+M2 


EBHF 




0.085 


0.2 


-33.99 


-18.05 


15.94 


-30.44 


-14.18 


16.26 


-23.00 


-6.92 


16.08 


16.22 


0.4 


-42.98 


-10.77 


32.20 


-39.75 


-7.01 


32.74 


-32.60 


-0.16 


32.44 


32.45 


0.6 


-51.71 


-3.31 


48.40 


-49.26 


0.15 


49.41 


-43.24 


5.67 


48.91 


48.67 


0.8 


-61.32 


3.68 


65.00 


-59.63 


6.58 


66.21 


-54.85 


10.49 


65.34 


64.90 


0.170 


0.2 


-45.94 


-23.19 


22.75 


-40.45 


-17.08 


23.37 


-31.01 


-8.24 


22.77 


23.00 


0.4 


-58.08 


-12.65 


45.43 


-53.46 


-6.72 


46.74 


-44.30 


1.25 


45.55 


46.00 


0.6 


-71.73 


-2.77 


68.96 


-68.10 


2.72 


70.82 


-59.60 


9.39 


68.99 


69.00 


0.8 


-86.22 


6.91 


93.13 


-83.99 


11.63 


95.62 


-76.42 


16.73 


93.15 


92.00 


0.340 


0.2 


-47.53 


-16.11 


31.42 


-38.89 


-5.75 


33.14 


-24.04 


7.12 


31.16 


32.30 


0.4 


-64.84 


-1.79 


63.05 


-57.48 


8.88 


66.36 


-42.31 


20.22 


62.53 


64.59 


0.6 


-82.75 


12.21 


94.96 


-77.17 


22.85 


900.02 


-62.44 


32.35 


94.79 


96.89 


0.8 


-103.15 


25.48 


128.63 


-99.87 


35.52 


135.39 


-86.32 


42.95 


129.27 


129.18 



5 Conclusions 

In this paper we have reported the study of asymmetric nuclear matter within the Brueckner- 
Bethe-Goldstone approach. The isospin effect on the equation of state has been investigated 
by performing a set of calculations at the two hole-line level of the BEG expansion for the 
energy per particle B{p,(3). The Bethe-Goldstone equation has been solved with the Ar- 
gonne V14 interaction. The continuous choice has been adopted for the auxiliary potential 
since it makes the convergence of the hole- line expansion faster than the gap choice P7| . 
Ranging the asymmetry parameter from j3 = (symmetric nuclear matter) to /3 = 1 (pure 
neutron matter) it was possible to check that B{p,(3) exhibits a linear dependence on 
for baryonic densities as large as at least two times the saturation density. This result con- 
firms the empirical law introduced in the mass formula of atomic nuclei and also extends 
its validity up to the highest asymmetries. As a consequence, the entire isospin effect is 
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incorporated in the symmetry energy. The calculation of the symmetry energy in the BHF 
approximation shows a monotonic increase as a function of baryonic density. Its value cal- 
culated at the saturation density is about 28.7MeV, in agreement with the empirical one. 
The comparison with the variational prediction is made rather difficult due to the contra- 
dictory results still existing in this approach. An accurate determination of the symmetry 
energy is required for dynamical simulations of collisions between neutron-rich nuclei, where 
the collective observables including collective flows, balance energy and other quantities are 
expected to be sensitive to the isospin degree of freedom 0. Ground state correlations 
were included in the mass operator up to the four hole-line order contributions. Their effect 
on the single particle properties has been investigated. The first-order contribution to the 
mass operator displays a linear dependence on the asymmetry parameter confirming a long- 
standing analysis by Lane |3^. A new effect of the isospin degree of freedom appears when 



the ground-state correlations induced by the second-order contribution are introduced in the 
mass operator. That is a nonlinear effect due to the particle-hole excitations of, say, protons 
induced by the propagation of a neutron in the nuclear medium. This new feature affects 
the isospin dependence of single-particle properties such as mean field, effective mass and 
mean free path. Along with the symmetry energy the heavy ion collisions with asymmetric 
nuclei could also probe the isospin dependence of mean free path and effective mass, which 
play also an important role in the collision dynamics. 

The EBHF approximation for asymmetric matter results in a satisfactory fulfillment 
of the Hugenholtz-Van Hove theorem in all asymmetry range < /5 < 1. This property 
makes us more confident of the hole-line expansion of the mass operator for calculating the 
single-particle properties including the Fermi energy. We found that the neutron and proton 
chemical potentials are largely affected by contributions beyond the BHF approximation. 
This could have far reaching consequences for the physics of the neutron star crust. In fact, 
the proton chemical potential in asymmetric nuclear matter is a very important ingredient 
in locating the inner boundary of the neutron star crust. However the difference fi = — /ip, 
and consequently the proton fraction in /3-stable matter, is almost unchanged in the EBHF 
approximation with respect to the BHF approximation. 
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